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HALO CONTRACTION EFFECT IN HYDRODYNAMIC SIMULATIONS OF GALAXY FORMATION 
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ABSTRACT 

The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark 
matter distribution. While this effect is based on simple gravitational physics, the question of its validity in 
hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state- 
of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from 
dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups 
with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We 
show that in all our cases the inner dark matter density increases relative to the matching simulation without 
baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to 
system and cannot be reduced to a simple prescription. We present a revised analytical model that describes 
the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the 
response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included 
in modeling of the mass distribution of galaxies and galaxy clusters. 

Subject headings: cosmology: theory — dark matter: halos: structure — galaxies: formation — methods: 
numerical simulations 



1. ADIABATIC CONTRACTION IN THE CONTEXT OF 
HIERARCHICAL GALAXY FORMATION 

Dissipationless cosmological simulations predict that viri- 
alized halos of dark matter are described by an approxi- 
mately universal density profile dD ubinski & Carlbera ll991t 
iNavarro et al][T997ll2010TlMoore et al.lll998l) . Although non- 
baryonic dark matter exceeds normal baryonic matter by a 
factor of ridm/^b ~ 5 on average in the universe, the gravita- 
tional field in the central regions of galaxies can be dominated 
by stars and gas. In the hierarchical galaxy formation picture, 
cosmic gas dissipates its thermal energy, condenses towards 
the halo center, and forms stars. In the process, dark matter 
particles are pulled inward and increase their central density. 

The response of dark matter to baryonic infall has tradition- 
ally been calculated using the model of adi a batic contraction 
(E22enetaD fT%l IZeldovich et ail [l980t iBarnes & White! 
1984). The present form of the sta ndard adiabatic contrac- 
tion (SAC) model was intr oduced by Blumen thal et al.l d 19861) 
and lRvden & Gunnl (119871) . This model assumes that a spher- 
ically symmetric halo can be thought of as a sequence of 
concentric shells, made of particles on circular orbits, which 
homologously contract while conserving the angular momen- 
tum. With these assumptions, the final location of the shell 
enclosing mass Md m (r), which was at radius r before the con- 
traction, can be calculated by knowing only the initial (Mb,,) 
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and final (Mb./) baryon mass profile: 

[M dm (r)+M b ,,(r)] r = [M to (r) +M hJ (r f )] r f . (1) 

Since the baryon mass typically increases at the inner radii as 
a result of dissipation (Mb./(r) > Mb,,(r)), the final radius of 
the shell is smaller than the initial radius, r/ < r. That is, the 
halo contracts. 

The effect of contraction of the dark matter distribution is 
important for many galactic studies: for modeling the mass 
profiles of galaxies and clusters of galaxies, for studying star 
formation feedback on the galactic structure, for abundance 
matching of galaxies and dark matter halos, for predicting 
possible signatures of dark matter annihilation, and others. 
This importance has led to a lively debate in the recent litera- 
ture, both theoretical and observational, on the validity of the 
effect (we discuss it in more detail in Section [2]). Hierarchi- 
cal galaxy formation is considerably more complex than the 
simple picture of quiescent cooling in a static spherical halo. 
Every halo is assembled via a series of mergers of smaller 
units, with the cooling of gas and contraction of dark matter 
occurring separately in each progenitor. Some objects may 
undergo dissipationless merging after the gas is exhausted or 
the cooling time becomes too long. Any of these effects may 
invalidate the assumptions of the SAC model and possibly its 
prediction of halo contraction. 

The only way to validate halo contraction is by inves- 
tigating cosmological hydrodynamic simulations that self- 
consistently model as many of these complex processes as 
po ssible: gas d i ssipat ion, star formation, and its feedback. 
In iGnedin et al.l d2004), we tested the SAC model using a 
suite of such simulations of the formation of one Milky Way- 
sized galaxy and eight clusters of galaxies, performed with 
the Adaptive Refinement Tree (ART) code. The comparison 
of the matching pairs of simulations with and without gas dis- 
sipation showed that the halos always contracted. However, 
the effect was weaker than that predicted by the SAC model. 
We developed a modified adiabatic contraction (MAC) model 
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based on the modified invariant, M(r)r, where r is the orbit- 
averaged radius for particles currently located at radius r. 

[MtoCP) +M bii (P)] r = [M Am (r)+M b , f (r f )] r f . (2) 

Although this combination is not strictly conserved during 
galaxy formation, it is a simple correction that most accurately 
predicts the dark matter profiles in the simulations. Using the 
mass within the average radius approximately accounts for 
eccentricity of particle orbits in a cosmologically-assembled 
halo. Averaging over particle orbits in radial bins, we found a 
mean relation between f and r in the range 10~ 3 < rj r vu - < 1 : 



' vir \ ' vir 



(3) 



where r vu - is the virial radius. This power-law dependence re- 
flects typical energy and eccentricity distributions of particles 
in cold dark matter halos in the dissipationless simulations. 
The parameters A and w varied from halo to halo and from 
epoch to epoch. For the systems at hand we found the mean 
values A w 0.85 and w w 0.8, which we used as fiducial pa- 
rameters of the model. We showed that the MAC model pre- 
dicts the halo mass profile to 10%-20% accuracy. 

Since 2004, several numerical and observational studies 
have evaluated the accuracy of the MAC model. We dis- 
cuss th ese studies in Section|2] In particular, Gustafsso n et al.1 
(2006) found different amounts of contraction in their four 
simulated halos and suggested that different combinations of 
the model parameters may best fit individual profiles. Fol- 
lowing the idea that the MAC model can be directly tested, 
in this paper we assemble a large collection of hydrodynamic 
simulations, of different systems run with different codes, to 
study the systematics of the contraction effect. We consider 
the effect of baryons only on the spherically-averaged radial 
halo distribution, although the halo shape i s also affected (e.g ., 
iKazantzidis eTani2004tlAbadi et aLlHoToHZemp et al1l20lll) . 



2. DEBATE IN THE LITERATURE 
2.1. Theoretical Studies 

Early con trolled simulations confirmed the halo contrac- 
tion effect (ISellwood & M cGaughl l200l IChoi et al.1 120061: 
IColm et al.l l 2006l). 

Recently, Duffv et al. (2010) considered 67 galactic and 
galaxy group halos in three cosmological simulations with 
different feedback prescription with the Smooth Particle Hy- 
drodynamics (SPH) code Gadget. They found the contraction 
effect in all cases, but each with a differen t amount of contrac- 
tion. Similarly to Gustaf sson et al.l d2006l) . they explored the 
best-fit distribution of parameters (A , w) and found them to be 
in the range 0.2-0.4 at z = 0. These parameter values may 
h ave been systematic ally lowered by the details of the fitting. 

Aba di et al.l (1201 Oh studied 13 galaxy-sized halos simulated 
with the SPH code GASOLINE. The simulations included ra- 
diative cooling of cosmic gas above 10 4 K but no cooling be- 
low 10 4 K and no star formation. They found that in each case 
the halo contracted but by a smaller amount than predicted by 
either SAC or MAC model. However, the weaker contrac- 
tion may be an artifact of the assumed gas physics, because 
their Fi gure 8 shows that the earlier simulat ions by the same 
authors (lAbadi et al.ll2003HMeza et al.ll2003l) . which included 
star formation and feedback, agree much better with the MAC 
m odel. 

iPedrosa et al.1 j2009ll20Toh and lTissera etail d20Tob studied 
6 halos in high-resolution SPH simulations with Gadget, with 



different prescriptions for star formation and feedback. They 
found halo contraction in all cases but again weaker than pre- 
dicted by either SAC or MAC model. They suggest that the 
dynamical formation history may affect the amount of con- 
traction in each individual halo. 

In contrast to these studies , in a simulation of 3 galaxies 
with Gadget, Johans son et al.l (120091) found that the dark mat- 
ter mass contained with a fixed radius of 2 kpc increased at 
high redshift, reached a peak at z ~ 3, and then declined con- 
tinuously to z = 0. The central density decreased by up to 
a factor of two between z = 3 and z = 0, reaching the same 
value as in a matching dissipationless simulation and thus ef- 
fectively canceling the contraction. Johans son et al.l (120091) 
attribute the reduction of the dark matter density to gravita- 
tional heating by infalling dense stellar clumps. In their run 
most of the baryons in the halo were converted into stars at 
high redshift (stellar mass is 13% of the virial mass), creat- 
ing massive stellar satellites, which may be too effective at 
pushing out dark matter. The evolution of the central den- 
sity may also be affected by a spurious numerical effect due 
to two-body scattering of massive particles. This scattering 
would alter the profile on the two-body relaxation time, which 
at 2 kpc is about 9 Gyr in their highest-resolution run, and is 
shorter in other runs. 

The SPH simulation of a large galaxy by 
iRomano-Dfaz et aT] d2008l) with constrained, rather than 
cosmological, initial conditions yielded an even more un- 
expected result. While they found the steepening of the 
dark matter profile at high redshift, the profile subsequently 
flattened to a constant-density core within the inner 3 kpc 
by z = 0, reversing a contraction into an expansion. Such 
flattening is likely to be a numerical relaxation effect. Based 
on the information in their paper, we estimate that there are 
about 9000 particles within 3 kpc, which have a two-body 
relaxation time of about 14 Gyr. Thus the interaction of 
these particles can alter the true profile within the age of the 
universe. 

Two new GA SOLINE simulations, one for a dwarf galaxy 
( Governato et al.l20 1 Oh and one for a Milky Way-sized galaxy 
(Guedes et al. 201 1), reach different conclusions on halo con- 
traction. In the first case the central dark matter density is 
reduced, while in the second it is enhanced. Both simulations 
use the same strong blastwave feedback, which suppresses 
gas cooling around the star-forming region for several Myr 
and creates powerful outflows that remove a significant frac- 
tion of baryons from the halo (70% in the first run and 30% 
in the second run). The difference can be tr aced to a higher 
density threshold for star formation in the iGovernato et al.l 
(2010) simulation, which means the supernova energy re- 
lease is more concentrate d and creates rapid potential fluctua- 
tions near the center. As lPontzen & Governato! d201 lb empha- 
size, the decrease of the inner dark matter density is achieved 
by repeated fluctuations following bursts of star formation, 
because a sing le outflow event ca nnot cause a significant 
change (e.g.. iGn edin & Zhao 2002). However, note that the 
IGovernato et al.l d2010l) result depends directl y on the adopted 
param eters (their lower-threshold version in [Governato et al.l 
120071 showed halo contraction instead of expansion) and thus 
it is one of many possible outcomes. 

Overall, cosmological simulations performed by different 
authors with very different codes and physics input agree that 
the contraction effect is present, but at a weaker level than 
suggested by the SAC model and with a significant variation 
from system to system. 
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2.2. Observational Studies 

Observational studies have general ly provided evidence for 
halo contraction. Schulz et al. (2010) studied the mass distri- 
bution of a stacked sample of 75,000 elliptical galaxies from 
the Sloan Digital Sky Survey, combining weak lensing analy- 
sis in the outer regions of the halo with measurements of the 
stellar velocity dispersion in the inner regions. They found 
that halo contraction is required to explain the large and sig- 
nificant mass excess over the dissipationless NFW profile, in 
all luminosity bins. The MAC model gave excellent predic- 
tions to the data, while the SAC model overestimated the ex- 
cess. An alternative explanation of the observed mass excess, 
by changing the stellar mass-to-light ratio, would require dou- 
bl ing the stella r mass relative to the Rroupa IMF. 

Aug er et all (120101) studied the mass distribution of 53 
early-type galaxies from the Sloan Lens ACS Survey, using 
a combination of stellar velocity dispersion measurements, 
strong lensing, and weak lensing. They found that the data 
favor a heavy stellar IMF and that halo contraction predicted 
by the MAC model improves the fit relative to the SAC model 
o r to the case without halo con traction. 

Minor & Kaplinghat (2008) investigated the effect of halo 
contraction on the statistics of strong gravitational lensing. 
They find that the contraction effe ct enhances the to tal lens- 
ing probability (similarly found by iRozo et a l. 2008) and in- 
creases the ratio of double images relative to quad images and 
lensing cusps. In particular, without contraction naked cusp 
configurations would become dominant at angular separation 
as small as 2.5", which is inconsistent with the data. 

While earl y X-ray studies of the m ass distribution in galaxy 
cluste rs (e.g..lZa ppacosta et al. 2006) disfavored halo contrac- 
tion, iDemocles et al.1 (120101) studied the mass distribution in 
fossil galaxy groups (otherw ise known as "X-Ray Overlumi- 
nous Elliptical Galaxies", see lVikhlinin et alJ 19991) and found 
th at halo contraction slightly improves the model fit. 

Seig ar et ail ((2008) studied the mass distribution in the M3 1 
galaxy and found that halo contraction is required to fit the 
inner regions. Without contraction, they cannot find any set of 
halo parameters that would reproduce the peak of the optical 
rotation curve at 10 kpc from the center, even with a very large 
halo mass. In fact, the best fit requires the SAC model, while 
the fiducial MAC model would need an unreasonably large 
c oncentration parameter o f the dark matter halo. 

Napolitano et aT] (1201 ll) derived the mass profile of a giant 
elliptical galaxy M84 using radial velocity measurements of 
planetary nebulae and field stars. The inclusion of the halo 
contraction effect, as parametrized by the MAC model, has 
allowed to derive the best-fit parameters of the dark halo to 
be in excellent agreement with the predictions of cosmolog- 
ical simulations. It is the first time that such agreement was 
a chieved with t h e plane t ary ne bulae kinematics. 

Dutt on et al.1 (120071 1201 ll) constructed semi-analytical 
models that simultaneously fit the luminosity function of 
galaxies and the Tully-Fisher/Faber-Jackson scaling relations. 
They found a preference for halo contraction for early-type 
galaxies, but a need for the opposite (halo expansion) for late- 
type galaxies. This conclusion is subject to additional uncer- 
ta inty in the stellar mass-t o-light ratio or the stellar IMF. 

Ben son & Bowed (2010) investigated extensive observa- 
tional data (galaxy sizes, clustering, luminosity function, 
mass-metallicity relation, Tully-Fisher relation, cosmic star 
formation history) using an improved version of the semi- 
analytical code GALFORM and a flexible prescription for halo 



contraction using the MAC model. Their best-fit model re- 
quired strong contraction (A = 0.74, w = 0.92; see also their 
Figure 7). Overall, they found a range of MAC model param- 
eters (A = 0.74-0.96, w = 0.81-0.99) matching the multiple 
data constraints. 

3. WHAT CONTROLS THE AMOUNT OF 
CONTRACTION? 

iGnedin et al.1 (|2004) suggested that the eccentricity of dark 
matter particle orbits may be responsible for the weaker con- 
traction effect seen in the simulations, relative to the SAC 
model. During the process of baryo n condensation p article 
orbits can further deform, although Debattista e t alJ (12008b 
find that this effect is almost completely reversible if the 
baryons are subsequently removed (for example, by strong 
stellar feedback). In addition, during violent changes of the 
gravitational potential, such as those during galaxy mergers 
or interactions, the orbital structure may differ from that in a 
dynamical equilibrium. Then the orbital action variables are 
not strictly conserved and a first-principles calculation of the 
contraction effect becomes even more difficult. 

Our original power-law approximation for the average orbit 
radius f was derived as an average of several dissipationless 
simulations of galaxy-size d and cluster-sized halos . Hydrody- 
namic SPH simulations of Gustafss on et alJ (j2006j) confirmed 
that Equation (O is an excellent fit in the inner parts of halos. 
They fit each of their four simulated halos with an indepen- 
dent set of the parameters A and w , and found them to lie in the 
range A = 0.74-0.83, w = 0.69-0.81. The values of the pa- 
rameters in our MAC model are consistent wi th these ranges. 
However, the results of Gustafsson et al. (2006) indicated that 
each simulated system may have somewhat differing amounts 
of contraction. 

To better understand the orbital dependence, we investi- 
gated the r—r relation in our hydr odynamic run for " Clus- 
ter 6" (shown in Figures 1 and 2 in lGnedin et al.ll2004l) . We 
computed a proxy for the average orbital radius, r ~ (r), by 
integrating orbits in a static NFW potential normalized to the 
maximum circular velocity V max measured in the simulation. 
Note that while the normalization of the profile correctly takes 
into account the baryon condensation, the shape of the poten- 
tial may deviate from the NFW form in the dissipative runs. 
Figure|T|shows that while the dissipationless run is well char- 
acterized by our fiducial slope w « 0.8, the run with cooling 
and star formation shows a more complex relation, with de- 
partures from a single power law. The outer regions, between 
0.1 r v i, and r v ; r , have a steep slope w 1 (the virial radius is 
calculated at the overdensity of 180 with respect to the mean 
cosmic density). The inner regions have m uch flatter slope, 
w ~ .5, similar to the values reported by Gustafsso n et alJ 
(2006). In these inner regions, r < 0.03r v i r , the contribution of 
baryons to the mass profile becomes significant and the con- 
traction effect should be most important. Thus for the purpose 
of calibrating the contraction model, we can restrict the fit of 
a power-law slope just to the inner regions. 

Interestingly, in the run in which cooling and star formation 
were turned off at z < 2, the f— r relation in the outer halo is 
closer to that in the dissipationless run than in the full dissipa- 
tive run. In the inner halo, the slope is again flatter, w « 0.6. 
In this run the stellar fraction within the virial radius is half 
of that in the full run and is closer to the observed fraction in 
galaxy clusters. After the cooling was stopped, the gravita- 
tional potential in the inner halo probably did not evolve and 
allowed the particle orbits to come to dynamical equilibrium. 
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which after dividing both sides by M,(r)r/ becomes 



FIG. 1. — Average orbital radius vs. instantaneous radius of dark matter 
particles in a dissipationless and two dissipational runs of a cluster-sized halo 
at z = 0. This halo is described in the second line in TablefT]and in Section l4~T1 
The size of the smallest resolution element is A = 1.7 X 1(T 3 r v i r . In one sim- 
ulation the gas cooling is suppressed at z < 2. Dotted and dashed lines show 
the slopes in the inner halo, with w = 0.8 and w = 0.6, respectively. Long- 
dashed line shows the relation in the SAC model. Vertical dotted line marks 
new pivot radius ro in Equation|4] Vertical arrow illustrates an adjustment by 
factor of 1.6 required for new normalization Ay. 



Therefore, the contraction effect may depend not only on the 
final amount of baryon condensation but also on the duration 
of this process. 

Since the f— r relation evolves during the process of dissi- 
pation, a MAC model with fixed parameters A and w cannot 
accurately describe the contraction effect. Instead, following 
Gustafsson et al] (120061) . we can search for the best-fitting val- 
ues of the model parameters for each simulated halo and then 
analyze their distribution. 

Figure Q] also shows that if we are to vary w, the normal- 
ization A would have to vary correspondingly because Equa- 
tion (0) is anchored at r v i r , far from the region of interest in 
the inner halo. We can reduce this correlated variation of A 
by shifting the pivot of the r— r relation to some inner radius, 
rn. The correlation between w and A, which we derive for the 
simulations described below, is minimized for rn w 0.03 r wu . 
Therefore, we redefine the relation of the MAC model as 



(4) 



We use the subscript "0" to differentiat e this normalization 
parameter Ao from its original version in lGnedin et al.l (2004) 
defined by Equation (0). We fix the pivot at ro = 0.03 r v i r in all 
discussion that follows. Only the normalization Ao is affected 
by this shift of the pivot point from r v h- to ro, the slope w 
remains invariant. The SAC model is still characterized by 
Ao = w = 1 . 

3.1. The role of model parameters 

To understand the effect of varying the parameters A and 
w on the amount of contraction, consider the inner regions of 
a halo. The contraction factor y = rj/r obeys Equation (|2), 



y 



= l-fb + 



Mb(ry) 
Mi(r) ' 



(5) 



where fb is the baryon fraction within the virial radius of the 
halo of interest. The average cosmic baryon fraction is fb ~ 
0.17, but a given halo may have a different baryon fraction. 
Here M,(r) = M dm ,,(r) + M b ,,(r) is the total initial mass, the 
sum of dark matter and baryons. Typically, before contraction 
the dark matter profile shadows the total mass distribution, 
M dm , / (r)«(l-/ 6 )M 1 -(r). 

In the inner halo described by an NFW profile with the scale 
radius r s , the initial enclosed mass is M,(r) oc r 2 at r <C r s . In 
this region, we can define the average logarithmic slope of the 
final baryon density profile: p n oc r~ v ', which lies in th e range 
l<z/<3(V«2is typical; Koopmans et al. 2009). Then 
the final baryon mass profile Mb(r) oc r i ~ v , and Equation (f5]l 
becomes 

1 . . M b (r) 



~ = 1 

y 



f b + 



M,(r) 



(6) 



The enhancement factor of the dark matter m ass profile is 
given by Equation (All) of iGnedin et al.1 (120041) and is most 
easily evaluated at the contracted radius ry = ry(r): 



F M (ry) 



M dm . f (ry) _ M dmJ (r) 1 



Afdm,i(r)0 M dm .i(ry) y{rf 



(7) 



The last approximation is valid only in the inner region where 
M,(r) oc r 2 . For arbitrary v and w, Equation (O is transcenden- 
tal and must be solved numerically. However, a good second- 
order approximation is described in Appendix. 

To make concrete calculations using Equations © and (0, 
we must specify the ratio of the final baryon enclosed mass to 
the initial total enclosed mass. Let r e be the radius where 
the final baryon mass equals the initial dark matter mass: 
Mb(r e ) = M dm j(r e ). For the Milky Way galaxy, this radius 
is about 10 kpc, and thus a good choice is r e w 0.05 r v ; r . It 
then follows that M b (r)/M;(f) = (1 - f h )(f / 'r e ) l - v '. 

Figure [2] shows the mass enhancement factor using the nu- 
merical solutions of Equation © for two representative val- 
ues of v = 2 and v = 1 .5. In the inner halo, the stronger baryon 
dissipation (v = 2 vs. v = 1.5) leads to the stronger contrac- 
tion of dark matter. Larger value of w and smaller value of Ao 
indicate stronger contraction effect at small radii. However, 
there is a cross-over point for the lines of different w, and 
the trend is reversed at the intermediate radii (0.01 -0.1 r v ; r ). 
While both Ao and w determine the normalization of Fm, the 
parameter w determines also the increase of the amount of 
contraction with radius. 

4. COSMOLOGICAL HYDRODYNAMIC SIMULATIONS 

In this section we describe the simulations of galaxy forma- 
tion that we use to test the MAC model. We tried to collect 
as diverse samples of simulations and codes as possible and 
to analyze them at z = 0, so as to evaluate the applicability of 
the model to observations. We also consider a few interesting 
cases at z = 1 and z = 3. 

We fit both the original set of parameters A and w from 
Equation ((3) and the modified set Ao and w from Equation 
©. The latter set for each system is summarized in TableQ] 

To derive the best-fit parameters we minimize the difference 
between the enclosed mass in a radial bin in the dissipative 
simulation, M s - lm , and the corresponding mass predicted by the 
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FIG. 2. — Enhancement of the enclosed dark matter mass given by Equa- 
tion (7) for the inner halo. The top shaded region corresponds to the baryon 
density slope v = 2, the bottom region to v = 1.5. Each region is bounded by 
the lines of the contraction models with the MAC parameter w = 1 (upper) 
and w = 0.6 (lower). The middle line is for w = 0.8. The parameter Aq is set 
to 1.6 in all cases. The baryon mass profile is normalized to be equal to the 
initial dark matter mass at r e = 0.05 r vlr . 



MAC model, M moc j. For convenience, we define the following 
dimensionless difference in radial bin z": 



Mr 



Zi 



1. 



(8) 



The "correct" mass profile in the dissipative simulation is not 
known infinitely precisely, and therefore the model should 
match it only within certain fidelity, cr,. The best-fit values 
of the parameters are obtained by minimizing the x 2 function 



X 



(91 



We include two sources of uncertainty in <t,. The first is a 
simple Poisson counting error, which we take to be the square 

-1 12 

root of the number of particles in the bin, N- t . The other 
is systematic uncertainty in the spherical mass profile, which 
can arise from a number of effects such as a triaxiality of the 
halo, an ongoing merger event, or a small mismatch in the out- 
put timing between the dissipational and dissipationless sim- 
ulations. In such cases the deviation in some radial bins be- 
tween the predicted model profile from the simulated profile 
may be large, even if the bin contains thousands of particles. 
We model this uncertainty as a constant error in each bin, cr sys , 
that puts an upper limit on the value of x~ ■ This allows us to 
obtain controlled, even if only relative, estimates of the confi- 
dence intervals of the model parameters. The value of cr sys is 
set such that for the SAC model (Aq = w= 1) it would result in 
a specified value of \ 2 per degree of freedom (equal to N sys ): 



T 2 = 
sys — 



l 



l 



N sys N b -2 



(10) 



TABLE 1 

Simulations and MAC model parameters 



System 




M vir (M Q ) 


An 


w 


rms 


»/" 


rms" 


Ref 


Cluster 




3.4 x 10 14 


1.83 


1.33 


0.037 


1.29 


0.060 


G04 


Cluster 




3.7 X 10 14 


1.74 


0.51 


0.028 


0.50 


0.045 


G04 


Cluster 




3.8 x 10 14 


1.54 


0.74 


0.015 


0.76 


0.022 


G04 


Cluster 




2.1 x 10 14 


2.30 


1.48 


0.051 


1.57 


0.140 


G04 


Cluster 




1.4 x 10 14 


1.78 


0.71 


0.009 


0.79 


0.053 


G04 


Cluster 




2.0 x 10 14 


1.28 


0.55 


0.027 


0.37 


0.100 


G04 


Cluster 




1.9 x 10 14 


1.54 


1.03 


0.007 


1.02 


0.016 


G04 


Cluster 




1.2 x 10 14 


1.24 


0.87 


0.013 


0.61 


0.232 


G04 


Group 




8.6 x 10 13 


1.42 


0.67 


0.009 


0.55 


0.032 


N06 


Group 




4.0 x 10 13 


1.30 


1.06 


0.005 


0.81 


0.030 


N06 


Group 




7.1 x 10 13 


1.47 


0.91 


0.011 


0.82 


0.024 


N06 


Group 




1.1 x 10 14 


1.43 


1.11 


0.023 


0.98 


0.036 


N06 


Group 




5.3 X 10 13 


1.14 


1.22 


0.035 


0.84 


0.067 


N06 


Group 




7.3 X 10 13 


2.06 


0.91 


0.012 


1.24 


0.078 


N06 


Group 




9.8 X 10 13 


1.59 


0.85 


0.017 


0.84 


0.057 


N06 


Group 




7.6 X 10 13 


1.41 


0.69 


0.009 


0.55 


0.031 


N06 


Group 




6.1 x 10 13 


0.75 


1.28 


0.008 


0.59 


0.135 


N06 


Group 




2.1 x 10 13 


1.12 


1.29 


0.020 


0.90 


0.048 


N06 


Group 




1.9 x 10 13 


1.87 


0.62 


0.010 


0.77 


0.017 


N06 


Group 




1.2 x 10 14 


1.61 


0.86 


0.027 


0.86 


0.027 


N06 


Galaxy 




5.5 x 10 11 


2.65 


1.38 


0.009 


0.85 


0.107 


G10 


Galaxy 




7.5 x 10 11 


1.79 


1.20 


0.014 


1.07 


0.034 


G10 


Galaxy 




2.8 x 10" 


1.21 


0.67 


0.016 


0.91 


0.068 


G10 


Galaxy z = 


1 


8.4 x 10" 


2.07 


0.64 


0.021 


0.99 


0.081 


C09 


Galaxy z = 


1 


4.8 x 10 10 


2.92 


0.85 


0.023 


1.31 


0.266 


C09 


Galaxy z = 


3 


3.4 x 10 11 


1.32 


1.26 


0.103 


1.26 


0.144 


L08 



Note: " for fixed Aq 

Refer ences: G04 

= IGottl oeber et al] 
Levine et al. (2008). 



= 1.6. 

= IGnedin et all f2004l). N06 = INagail (2006), G10 
(2010), C09 = ICeverino & Klypin (2009), L08 = 



freedom is A^-2 when we are fitting two model parameters. 
Once we set cr sys , it is then kept fixed as we search for the 
best-fitting parameters A and w. 

We take a large enough value, N sys = 100, such that it would 
not alter fitting in the bins with fewer than 100 particles, while 
avoiding large deviations in the more populous bins. We have 
verified that choosing any value in the range 10 < N sys < 1000 
results in similar model parameters. The total error in bin i is 



^- 1 + „ 2 

Uj = 1" <7„.,„ . 

' Ni sys 



(11) 



where Nb is the number of bins, and the number of degrees of 



As an alternative to minimizing the \ 2 function, we have 
also tested a robust estimator ^ \Zi\/&u which is not sensitive 
to distant outliers. With both methods we obtained essentially 
the same best-fit parameters for our simulations. 

In order to make the most accurate fit in the region where 
the contraction effect is important, we restrict the radial bins 
included in the fitting to r < 0.1 r v ; r . We allow the parameters 
to vary in the range 0<Ao<3,0<w<2. Even larger val- 
ues would yield models degenerate with those in the chosen 
range, as it will become clear later in Figure [6] 

4. 1 . Clusters and groups of galaxies 

The first sample consi sts of the simulatio ns of 8 galaxy 
cluster halos described in Gnedin et al. (2004) and additional 
12 galaxy group halos by Nagai (2006) with the same setup. 
These are high-resolution cosmological simulations in the 
ACDM model (fl m = 0.3, fl A = 0.7, Q b = 0.043, h = 0.7 
cr 8 = 0.9) performed with the ART code dKravtsovl Il999t 
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FIG. 3. — Contractio n of the dark matter profile in a simulated group of 
galaxies at z = 0, from Nagai (2006). Solid lines show the enclosed dark 
matter mass profile, in the non-radiative (ad) run and star formation (sf) run. 
Dotted lines show the corresponding baryon mass profiles. Black solid line 
is the best fit of the MAC model, with parameters Aq = 1.61, w = 0.86. Top 
panel shows the mass residuals for the MAC model with freely adjustable Aq 
(solid) and the MAC model with fixed Aq = 1.6 (dashed). In this plot the two 
lines almost coincide. 



iKravtsov. Klvpin. & Hoff man 2002). The simulations have 
a peak spatial resolution A = 3.5 kpc and dark matter particle 
mass of 3.9 x 10 8 Mq. The virial mass of the systems ranges 
from 2 x 10 13 M Q to 4 x 10 14 M Q . Star formation is imple- 
mented using the standard Kennicutt's law and is allowed to 
proceed in regions with temperature T < 10 4 K and gas den- 
sity n g > 0.1 cm" 3 . We truncate the inner profiles at 4A to 
ensure that the gravitational dynamics is calculated correctly 
in the studied region. 

Figure [3] shows the mass profiles for one of the groups. The 
dark matter mass is significantly enhanced in the star forma- 
tion run relative to the non-radiative run, by a factor of 4 at 
the innermost resolved radius. The baryons strongly domi- 
nate the total mass at that point. The MAC model provides 
an excellent fit to the contracted dark matter profile, with the 
parameters {Aq = 1.61, w = 0.86) close to the fiducial values. 
The maximum deviation of the mass profile predicted by the 
MAC model is 6%, and the rms deviation over all bins at radii 
r < 0.1r v j r is 3%. We similarly analyzed the other eleven 
groups and present their best-fit parameters in the discussion 
of Figure [6] 

4.2. Individual Galaxies 

We consider the simulation of three Milky Way-sized 
galaxies by the CLUES project (http://ww w. clues-proj ect. org ; 
Gottloeber et al] l2010t lKnebe et al. 120101) . The simulation is 
run using the SPH code Gadget-2. This code includes stan- 
dard radiative cooling, star formation, and supernova feed- 
back. The force softening length e = 0.14 kpc. The halos were 
selected from a large box and resimulated with the effective 
mass resolution of 4096 3 dark matter particles. In the highest- 
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r (kpc) 

FIG. 4. — - Contraction of the da rk matter profile in a simulated galaxy at 
z = 0, from Gottloeber et al. 12010). Line notation is the same as in FigurefJ] 
Black solid line is the best-fit MAC model with A = 1.79, w = 1.2. Dashed 
line in top panel shows the MAC model prediction with fixed Aq = 1.6, and 
best-fitting w = 1.07. 



resolution halos the particle mass is 3.5 x 10 5 M Q . The virial 
masses of the three halos at z = are (3 - 8) x 10 11 Mq. The 
inner truncation radius is set by the condition that the local 
two-body relaxation time exceeds the age of the universe. 

Figure 2] shows the profile of the most massive of the three 
galaxies. The dark matter mass is enhanced by an order of 
magnitude at the innermost radius. The MAC model with pa- 
rameters (Aq = 1.79, w = 1.2) predicts the dark matter profile 
to better than 4% accuracy in any bin, with the rms deviation 
of only 1.4%. 

We consider also the simulations of a Milky Way-sized 
galaxy and a dwarf galaxy at z = 1 by ICeverino & K lvpin 
(2009). These simulations are run with the ART cod e with 
a very different prescription for stellar feedback than in lNagail 
(2006). The large galaxy mass is 8 x 1O U M , the dwarf 
galaxy mass is 5 x 1O 1O M , both at z = 1. The dark matter 
particle mass is 7.5 x 10 5 M Q and the peak spatial resolution is 
100 comoving pc for the larger galaxy. For the smaller galaxy, 
the dark matter particle mass is 9.4 x 10 4 M Q and the peak 
resolution is 50 comoving pc. Compared to the non-radiative 
runs, the dark matter mass is enhanced by a factor of 8 for 
the larger galaxy and by a factor of 5 for the smaller galaxy, 
at the innermost radius. The MAC model (with parameters 
A = 2.07, w = 0.64 and Aq = 2.92, w = 0.85, respectively) pre- 
dicts the dark matter profile to better than 9% accuracy, with 
the rms deviation of about 2%. 

4.3. Galaxy Center 

Finally , we consider the res imulation of the galaxy run re- 
ported in iGnedin et alj d2004l) that zooms into th e innermost 
region of the galaxy at z = 3 dLevine et alj|2008l) . This sim- 
ulation follows the early evolution of a galaxy that becomes 
a Milky Way-sized object at z = 0. The DM particle mass is 
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FIG. 5. — Contraction of the dark matter profile in a simulated galaxy at 
7, = 3, from Levine et al. 12008). Line notation is the same as in Figure [3] 
Black solid line is the best-fit MAC model with A = 1.32, w = 1.26. Vertical 
bars on the MAC model in the top panel indicate the Poisson uncertainty of 
the mass profile derived in the simulation. Dashed line in top panel shows the 
MAC model prediction with fixed Ag = 1.6, and best-fitting w = 1.26. 




FIG. 6. — Best-fitting parameters of the original MAC model (EquationfJ) 
for all simulations discuss ed in this paper. Aste risk marks the fiducial param- 
eters of the MAC model in Gnedin et al. 12004). Solid lines show the relation 
between A and w that gives the same amount of contraction (enhancement 
of dark matter mass) at r = 0.005 r v j r , for the baryon profile with v — 2 nor- 
malized to equal the initial dark matter mass at r e — 0.05 r v ; r . The top line 
gives the same amount of contraction as the SAC model. The other two lines 
correspond to 50% and 30% of that amount. 



1.3 x 1O 6 M and the peak force resolution at z = 3 is 0.064 
kpc for the gas and 0. 1 kpc for the dark matter, a very small 
scale for cosmological simulations. We truncate the inner pro- 
file such that the innermost bin contains at least 200 dark mat- 
ter particles. 

Figure|5]shows that the MAC model is able to describe even 
this case, with the rms deviation of 10%. This case is extreme 
because the baryons dominate the dark matter by two orders 
of magnitude at the innermost radius, and the dark matter 
mass is enhanced by a factor of 300 relative to the extrapo- 
lation of the dissipationless profile. 

We also note that the stellar profile is contracted similarly 
to the dark mater profile, because gas accretion is faster than 
star formation. 

5. DISTRIBUTION OF MODEL PARAMETERS 

All of the simulations considered here indicate some de- 
gree of enhancement of the dark matter profile. Not a single 
case indicates halo expansion rather than contraction. Fig- 
ure [6] combines the resulting constraints on the parameters A 
and w of Equation (01. The models do not fill all the available 
parameter space, but instead concentrate in a fairly narrow re- 
gion in which A and w are strongly correlated. The original 
MAC model suggested by iGnedin et al.l (120041) falls right in 
the middle of the new distribution. 

It is interesting to determine which combination of the pa- 
rameters A and w yields the same amount of contraction. 
Given the radial dependence of the mass enhancement fac- 
tor Fm (Equation 17), the solution to this problem varies with 
radius. However, we can remove most of the radial depen- 
dence by defining the enhancement factor relative to the SAC 
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FIG. 7. — Best-fitting parameters of the revised MAC model with rg ■■ 
0.03 r v j r (Equationf4). Symbols and lines are as in Figure|6] 



model: 



= F M (r\A,w) 
■ tM ~ FmO|1,1) 



(12) 



and evaluating it at some inner radius where the linear ap- 
proximation for the contraction factor y(r) is valid. We take 
r = 0.005 r v h, which corresponds to about 1 kpc for the Milky 
Way galaxy. The exact value of r affects the resulting value 
of parameter w (for a given A) only logarithmically, as long as 
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FIG. 8. — Comparison of the best-fitting parameters w obtained with vari- 
able Aq with those obtained for fixed Aq = 1.6. Sizes of the errorbars are 
relative. 
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FIG. 9. — Best-fitting parameter w (for fixed Aq = 1.6) vs. halo mass for all 
systems in our compilation. 



r r e , which we take again to be r e = 0.05 r vu -. 

Lines in Figure|6]show the relation between A and w corre- 
sponding to three values of /m- All simulations but three fall 
below the level of contraction predicted by the SAC model. At 
the same time, no simulation falls below the level of fu = 0.3. 
Therefore, the MAC model is well constrained to be able to 
reliably predict the amount of dark matter in the inner regions 
of galaxies and clusters. 

The isocontours of constant /m become even more horizon- 
tal at smaller r. This suggests that parameter w may be more 
important than A in describing the amount of contraction. It 
is also desirable to describe the strength of the contraction 
effect by only one parameter instead of two. The first step in 
this direction is to eliminate or reduce the apparent correlation 
between A and w. To this aim, we calculated the best-fitting 
parameters of the revised model (Equation|4]i and determined 
the pivot radius ro that minimizes the correlation. It is the 
value ro = 0.03 r vu - quoted above. Figure [7] shows the new 
distribution of the resulting best-fit parameters Aq and w. The 
values of w are essentially unchanged, but the values of Aq are 
more concentrated than the distribution of A in Figure|6] The 
residual scatter of Aq reflects intrinsic variation of the strength 
of halo contraction among different systems. The isocontours 
of constant /m also correspondingly change shape. 

The second step in simplifying the model prescription is 
fixing the parameter Aq. We take the average value Aq « 1.6 
and redo all model fits allowing only for the variation of w. 
The best-fit values are listed in Table Q] Obviously, a one- 
parameter fit is less accurate than the two-parameter fit, but 
the rms error of mass is still typically below 10%. The one- 
parameter fits also do not introduce systematic shifts in the 
derived values of w, as shown in Figure [8] Thus, w can serve 
as a convenient measure of the strength of halo contraction. 

In practical application of the contraction model to observa- 
tions or dissipationless simulations we wish to know the value 
of w most appropriate to a given system. We considered sev- 
eral properties of the simulated halos but, unfortunately, we 
were unable to find significant correlation with w. For exam- 
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FIG. 10. — Best-fitting parameter w (for fixed Aq = 1.6) vs. the bulge frac- 
tion, defined as the fraction of baryon mass contained within 0.01 r vu -, 



pie, Figure [9] shows that w is effectively independent of halo 
mass. Only the lower envelope of the distribution decreases 
with M vir . 

We have looked for other potential correlations: with the 
ratio of final baryon mass to initial dark matter mass at the 
innermost resolved radius and at a fixed radius of 0.01 r vu -, 
with the ratio of final dark matter to baryon mass at 0.01 r v j r , 
and with the bulge fraction of galaxies. The latter is defined as 
the fraction of baryon mass contained within 0.01 r wn ■: /bulge = 
M*(0.01 r v i r )/M*. Figure [TOl shows the scatter plot of w with 
the bulge fraction, lacking any significant correlation. 

6. IS HALO CONTRACTION REAL? 
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FIG. 11. — Model-independent strength of the contraction effect. Solid 
lines are for all systems in our c ompilation. Dashed blue line is the relation 
suggested by Abadi et al. (2010), dotted lines are the relations suggested by 
Pedrosa et al. 12010), dashed red line is a relation of the type suggested by 
Dutton et al. (201 1) but with n = 0.5 instead of n = -0.5 (which would go off 
the scale of the plot). 

The strength of the contraction effect can be assessed in a 
model-independent way, as suggested by I Abadi et al.1 (120101) . 
They expressed the ratio of the contracted to initial radius of 
spherical shells as a function of the ratio of initial to final total 
mass. In the SAC model this relation is expected to be unity, 
r $ I r = Mi(r) I M f(r f). For a weake r contr action effect, we ex- 
pect r f /r > Mi/Mf. lAbadi et all (120101) found that in their 
simulations the dark matter sh ells contract on the average as 
r f /r = 1 + 0.3 [(Mi/Mj) 2 - 1 ] . iPedrosa et all (120101) did a simi- 
lar analysis for their simulations but found different relations: 
r//r = 1 +0.14[(M,/M/) 4 - 1.25] for early-type galaxies, and 
r f /r= l+0.15[(M,/M/) 4 -1.4] for late-type galaxies. 

FigureQT]shows the contraction factor for all of our simula- 
tions. The effect is clearly more significant than suggested by 
the above relations. In several cases, the contraction is even 
stronger at the innermost radii than that predicted by the SAC 
model, that is rj/r < Mi/Mf. It can also be seen in Figure|7] 
Overall, there is no single relation for the contraction factor 
rf/r. The scatter among different systems is intrinsic, similar 
to the scatter of mod el parameters A and w. 

iDutton et al.l (1201 lb constructed semi-analytical models of 
galaxy populations using, in addition to the MAC model, an 
analytic prescription for either halo contraction or expansion 
of the form rj/r= (Mi/Mf)", with n as a free parameter. Their 
preferred case is halo expansion with n = -0.5. None of our 
systems shows evidence for such expansion. If we were to 
use this parametrization to envelop the distribution of our sim- 
ulations, then most of our results would lie within the range 
0.6 < n < 1 .2. However, such parametrization has no physical 
basis and we do not recommend it. 

Is the halo contraction effect real? Much of the apparent 
controversy in the literature on the validity of the contrac- 
tion effect is due to the strict application of adiabatic invari- 
ants. In fact, most of the hydrodynamic cosmological sim- 
ulations are in agreement that the effect i s real and, at the 
same time, weaker than that expected in the Blumenthal et al. 



(1986) model. To avoid future controversy, we propose to 
abandon the term "adiabatic contraction" (reserving it only 
for the SAC model for historical reasons) and instead use the 
term "halo contraction". 

Note that all individual effects that were proposed to "re- 
verse" contraction (such as the rapid supernova winds, cold 
accretion, galactic bars, inspiraling of dense baryonic clumps 
by dynamical friction, etc.) are already included in self- 
consistent cosmological simulations and differ only in the 
specific implementation (numerical resolution, star formation 
prescription, feedback model). Isolated investigations of these 
effects therefore do not invalidate the conclusions we derive 
from the ensemble of the simulations presented here. One 
should not also attribute halo contraction only to dissipative 
galaxy formation an d contrast it with dissipat ionless accre- 
tion of satellites (e.g.. lLackner & Ostrik erl2010l) because both 
processes are taking place simultaneously. Hydrodynamic 
simulations are steadily improving in accuracy and will con- 
tinue to include and evaluate new effects, such as the re- 
peated potential fluctuations follo wing bursts of star forma- 
tion dPontzen & Gov ernato 2011), mechanical and radiative 
feedback of active galactic nuclei, etc. Our ability to model 
halo contraction will continue to evolve along with our overall 
understanding of galaxy formation. 

The MAC model does not presume a specific amount of 
contraction. It is a method for calculating the dynamical re- 
sponse of the dark matter halo to a given accumulated mass 
of baryons. The model is based on simple underlying physics 
and not just fitting the results of a particular simulation. The 
halo response depends directly on the change of the baryon 
profile relative to the dissipationless formation. As such, the 
response may be a halo expansion if more baryons are re- 
moved fro m the galaxy center than we r e the re initially. For 
example, Sommer-Larsen & Limo usin! (120101) gradually re- 
duced the stellar mass in their simulated galaxy cluster be- 
tween z = 2 and z = 0, and found a much reduced contraction 
effect at the end, as would be expected if the removed stars 
did not form at all. 

How should halo contraction be included in modeling of 
observed systems or in theoretical semi-analytical models 
based on dissipationless simulations? Direc t parametrization 
of the form suggested by Aba di et al.1 (120101) can only be used 
when the final distribution of dark matter is known, that is, 
only for comparing results of hydrodynamic simulations with 
each other. When the distribution of dark matter needs to 
be predicted, the MAC model provides a reasonably accurate 
method and requires knowing only the final baryon profile. 
Unfortunately, the values of the model parameters have irre- 
ducible scatter and do not appear to correlate with any obvious 
property of the system of interest. Thus the model prediction 
carries some irreducible uncertainty. In practical application 
it suffices to account for this uncertainty by varying only one 
of the two parameters. We advocate the choice of Aq and w 
(EquationHJi instead of A and w (Equation[3]l because the new 
parameters are effectively uncorrelated. For assessing the ef- 
fect of halo contraction on a given system, we suggest fixing 
Ao =1.6 and varying w in the range 0.6- 1.3. This prescrip- 
tion covers most of the parameter range seen in our simula- 
tions (Figure |7]l while providing the estimate of the dark mat- 
ter mass profile with rms accuracy of about 10% (TableHJ. 

What is the origin of the intrinsic scatter of the parameters 
Aq and w, or of the contraction factor rf/rl We can identify a 
number of processes in galaxy formation that could produce 
such scatter. For example, baryons in different galaxies con- 
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dense at different epochs and have different angular momen- 
tum profiles, both of which affect the compactness of the final 
baryon distribution. In addition, galaxies experience different 
amounts of merger activity, which can change the morphol- 
ogy of baryon distribution and can gravitationally heat dark 
matter. Given these and other factors, it would be surprising 
if the intrinsic scatter did not exist. It is not known at the mo- 
ment which of these effects are dominant. It would be useful 
to explore the origin of the intrinsic scatter in future studies. 

7. CONCLUSIONS 

We have evaluated the halo contraction effect using a large 
compilation of cosmological hydrodynamic simulations per- 
formed by different groups with different codes. In all the 
cases we considered, we find an increase of the dark matter 
density in inner regions of galaxies and clusters, relative to 
the matching dissipationless simulation. The halo contraction 
effect is real and must be included in the modeling of obser- 
vations and in the semi-analytical theoretical modeling. 

The contraction effect is weaker than predic t ed b y the adi- 
abatic contraction model of Blumentha l et al.l (11986b . How- 
ever, depending on the system and the final baryon distribu- 
tion, the inner dark matter density is still enhanced typically 
by a factor of several, and in extreme cases by two orders of 
magnitude. 

The revised MAC model offers a convenient and accurate 



way to estimate the effect of halo contraction. The distri- 
bution of model parameters cannot be reduced to a single 
number, but their range is well constrained: 1 < Aq < 2.2, 
0.6 < w < 1.3. We find that fixing the value of Ao = 1 .6 does 
not significantly degrade the accuracy of the predicted mass 
profile, relative to the two-parameter fit. We suggest varying 
w in the range 0.6- 1 .3 in order to bracket a possible response 
of the dark matter halo in a given system of interest. 

The revised MAC model is encoded in the soft- 
ware package Contra, available for download at 
http://www.astro.lsa.umich.edu/^ognedin/contra 
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APPENDIX 

ANALYTICAL APPROXIMATION FOR THE CONTRACTION MODEL 

An approximate analytical solution of Equation (0 for the contraction factory = rj/r can be obtained as follows. At r <C r e , 
the ratio M b (r)/Mj(r) 3> 1 and the term (1 — /&) can be neglected. The first-order approximation is therefore 



yo : 



M b (r) 



Mi{r) 



(Al) 



Then we can express the correct solution as y = yo(l + S), with 5 < 1 . Substituting this into Equation © and expanding the power 
term of (1+8), we obtain 

d-/^o 

l+w(3-v) + (l-f b )y 

The value of 6 is of the order yo, which means that this approximation is valid where yo < 1, or r < r e . 

To evaluate the mass enhancement factor F^(r) at a specified radius r, we must express the initial radius r as a function of 
the contracted radius ry(r). Considering only the first-order approximation y w yo, we find a power-law solution ry oc r a , where 
a = (1 +2w)/(l + w(3-v)). Then F M (ry) oc (ry) °" 2 . Retaining all the coefficients, we have 



F M (ry) '■ 



~M b (?y 


2 

1+2^ 


Mi(f) 





(!-/*)(- 



(A3) 



This approximation is valid at r <C r e . 
Using the same first-order approximation, we can derive the inner logarithmic slope of the contracted dark matter profile, 

pdm(r) oc r" 7 : 

7 = 7T2 ^ (A4) 

This slope was already derived as Equation (A12) in Gne din et al.l (|2004). This equation links the contracted dark matter slope 
with the slope of the baryon profile. It provides an accurate description of the dark matter profile in the hydrodynamic simulations 
described above, within the errors of calculation of w and v. 
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